Strong Nernst-Ettingshausen effect in folded graphene 
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We study electronic transport in graphene under the influence of a transversal magnetic field 
B(r) — B(x)e z with the asymptotics B(x — > ±oo) = ±Bo, which could be realized via a folded 
graphene sheet in a constant magnetic field, for example. By solving the effective Dirac equation, 
we find robust modes with a finite energy gap which propagate along the fold - where particles 
and holes move in opposite directions. Exciting these particle-hole pairs with incident photons 
would then generate a nearly perfect charge separation and thus a strong magneto-thermoelectric 
(Nernst-Ettingshausen) or magneto-photoelectric effect - even at room temperature. 

PACS numbers: 72.80. Vp, 78.67.Wj, 85.80.Fi. 



Introduction The Nernst-Ettingshausen effect [l[ de- 
scribes the generation of an electric current (or voltage) 
by a temperature gradient in the presence of a magnetic 
field. Such thermoelectric effects facilitate the direct con- 
version of thermal into electric energy and thus are of 
general interest. Obviously, the £ (charge), *}3 (parity), 
and T (time reversal) symmetries must be broken for 
such an effect to occur. One way to achieve this is a 
magnetic field in a suitable geometry: trajectories of op- 
posite charge carriers are bent to antipodal directions. 
However, the mean free path in usual materials is too 
short to generate an efficient charge separation in that 
way - at least at room temperature. For example, the 
classical cyclotron radius r = m e v / (q e B) of a free elec- 
tron at room temperature in a magnetic field B of one 
Tesla r = C(/xm) is much larger than the typical mean 
free path (in the nanometer range). Thus, the Nernst- 
Ettingshausen effect is strongly suppressed by multiple 
scattering events and dissipation etc. 

This motivates the study of graphene since this 

system offers a comparably long mean free path and 
a large electron mobility, a linear (pseudo-relativistic) 
dispersion relation at low energies (i.e., near the Dirac 
points), and a very large Fermi velocity vp ~ 10 6 m/s Q, 
see also @, [H ■ In this case, the pseudo-relativistic cy- 
clotron radius at room temperature in a magnetic field 
of one Tesla is much smaller (some tens of nanometers). 
In this regime, quantum effects should be taken into ac- 
count - even at room temperature jij. 

In the following, we consider folded graphene in a 
transversal magnetic field, see Fig. [TJ In principle, the 
folding of graphene has already been realized experimen- 
tally, see, e.g., 0, This set-up is advantageous since 
we avoid real edges in graphene which are typically not 
perfect and contain cracks or other defects which might 
induce scattering, coupling to vibrational degrees of free- 
dom, or further unwanted effects. Form a theoretical 
point of view, these edges can only be described in ideal- 
ized cases, e.g., via effective boundary conditions which 
then depend on the concrete realization (e.g., zigzag or 
armchair structure 11-13^. 
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FIG. 1. Sketch of the considered set-up. 



Eigen-modes We consider length scales (e.g., curva- 
ture radius of fold) far above the lattice spacing of 
graphene « 0.25 nm and energies of 1 eV or below. In 
this limit, we may describe the low-energy behavior by an 
effective Dirac equation in 2+1 dimensions (H = q e = 1) 
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with x^ 1 = [vpt, x, y\, where vp « 10 6 m/s is the Fermi 
velocity [HI]. The Dirac matrices 7^ = [a z , io~ v , — ia x ] 
acting on $ = [ipi , ^2] are related to the Pauli matrices 
a x ' v ' z . In the Landau gauge, the vector potential = 
[0, 0, A(x)] generates the magnetic field B(x) = d x A(x) 
with the asymptotics B(x — > ±00) = ±Bq. 

In view of the translation symmetry in t and y, we can 
make the separation ansatz for the modes 



(2) 



x, y) = exp {-iEt + iky} ^ E ' k (x) , 

arriving at the two coupled equations 

iv F [d x + k + A(x)]^' k (x) = Eyb E - k (x) 

iv F [d x - k ~ A(x)]i/jf' k (x) = E^f k (x) . (3) 

Hence, we can choose ip^' (x) to be real, for example, 
while ifj2' k ( x ) i s imaginary. We observe a particle-hole 
symmetry since replacing E — > —E and ip2' k ~^ ~' l l ) 2' k 
yields a new solution ^~ E - k = a z ^ E - k = (^ E ^ k )*. 

The two first-order equations ([3]) can be combined into 
one second-order equation 



vp[k + A{x) + B x ][k + A(x) - d x ]ipl 
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and analogously for if)^' with d x <-> —d x . This equa- 
tion can be cast into the form of a one-dimensional 
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Schrodinger equation Hkif>i' k = E 2 ifjf' k with the Hamil- 
tonian Hk — Vp(—d 2 + Vk) containing the effective po- 
tential V fc = [k + A(x)] 2 + A'(x). Since this Hamilto- 
nian is self-adjoint Hk = H k and the potential Vk has 
the asymptotics Vk(x — ¥ ±00) = 00, we get a complete 
set of discrete, orthonormal, and localized (in x) eigen- 
functions ipf' k (x) for every value of k. These modes are 
non-degenerate for each k, i.e., the energy bands E(k) do 
not cross 
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Due to the particle-hole symmetry, each 
of these eigen- functions ipi' k {x) corresponds to a pair of 
modes ty ±E ' k (x) of the original problem ([3]) with oppo- 
site energies. Furthermore, with T>k = k + A(x) — d x , 
we may write Hk = v \T^kP k which shows that Hk is 
non- negative (and thus E is real). In addition, Hk can- 
not have a zero eigen-value E = since the correspond- 
ing ij}i~°' k (x) must satisfy T>kipi~°' k = 0, which gives 
?pf~ ' k (x) oc cxpjfcx + J dx A(x)} and analogously for 
il>2~°' k (x). Due to the asymptotics B(x — > ±00) = ±B a 
and thus A(x — > ±00) ~ -Bo Mi this solution is not nor- 
malizable and thus Hk is strictly positive for any k. Ergo, 
the modes do always have a finite energy gap E ^ 0. 
Current The current density of the modes reads 



^E,k 



V F ^E,kl^E,k = V^ E>k f^ E ,k ■ (5) 

The zeroth component j — v-pp is simply given by the 
density p = |V>f'T + |V>2 | 2 - As one would expect, 
j x vanishes identically since ^1 ' {%) is real and ip 2 ' ( x ) 
imaginary, cf. Eq. (|3|). Using the same argument, the 
current density in y-direction simplifies to 
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From the triangle inequality (2|a&| < \a 2 \ + \b 2 \), we may 
infer \j v \ < vpp, i.e., the speed of the associated charge 
carriers is at most the Fermi velocity wp (as expected). 
The total current in y-dircction can be obtained by 



J y 



dx j l 



2u| 
' E 



dx^' k [k + A{x)]^' k , (7) 



where we have used v-p'Dk'4 , i' = iEip^'" from Eq. (|3|). 
For the lowest E 2 modes (for a given fc), i.e., the upper- 
most negative mode and the lower-most positive mode, 
the wave- function ip^' k (x) corresponds to the ground 
state of Hk and hence it is non-zero for all x (node theo- 
rem (l5| ) . Since one can repeat the same line of argument 
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for tps' ( x )> the integrand j v = —2ivFipi'^ip2^ 1S non- 
zero for all x and hence the current J v is finite. But other 
modes could have J y = at some fc-value. However, for 
large enough k > — A mul = —imii{A(x)}, the integrand 
in the above equation ipf' k [k + A(x)]ipf' k is positive for 
all x and thus the current has a finite value. 

Furthermore, the current J v is related to the slope 
dE/dk of the dispersion relation, i.e., the group velocity: 
Writing Eq. © as H E .k \*E,k) = E \V E ,k), we find 
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where we have used the normalization e ,k\^ E ,k) = 1- 
Together with Eq. ([7]) we find that particles with E > 
and holes with E < have the opposite current (and 
group velocity), i.e., all particles (with k > — A m i n ) move 
to the right and all holes move to the left. In this way, 
one obtains a (nearly) perfect charge separation. 

Asymptotics It is illustrative to study the two limit- 
ing cases k — > ±00. For large and positive k, the po- 
tential Vk can be approximated by Vk ~ k 2 + 2kA(x). 
Thus, to lowest order in fc, we obtain E w iupfc, i-e., 
these modes propagate with a speed close to the Fermi 
velocity. Going to the next order in fc, we may expand 
A{x) around its minimum at xq where the magnetic field 
B(x ) = vanishes A(x) sa A miu + A"(x )(x ~ x ) 2 /2 
and obtain harmonic oscillator eigen-functions centered 
at xq [assuming that A"{xq) 7^ 0]. Since the stiffness of 
the potential behaves as kB\xo), the modes are strongly 
localized around xq for large k and basically propagate 
along the xo-line where the magnetic field vanishes. For 
fixed and large fc, these modes have equidistant values of 
E where the distance scales with y/ B'(xq). 

For large and negative fc- values, the minima of the po- 
tential Vfe are given by A(x±) + k = and thus the 
modes are localized at large and nearly opposite values 
of x± ~ ±\k/B \ due to A(x — > ±00) ~ -BoM- I n this 
regime, A{x) is approximately linear and thus we recover 
the harmonic oscillator eigen-functions corresponding to 
the usual (pseudo-relativistic) Landau levels in a con- 
stant magnetic field Note, however, that the eigen- 
functions ipf' k (x) are linear superpositions of the Lan- 
dau levels centered at x + and x_ with the same energy 
E. In this limit, the eigen-energies E do not depend on 
fc anymore = ±i>p \J2B$n with n £ N) and thus the 
current J v also vanishes. Hence these modes are not so 
interesting for our purpose. 

Matrix elements Now we are in the position to study 
the excitation of particle-hole pairs by incident photons 
(in the infra-red or optical regime). In second quantiza- 
tion, the interaction Hamiltonian reads 



dxdy ^^A^ . 



(9) 



where the photon field operator A^ contains the cre- 
ation and annihilation operators k and d w .fc jCr for 
frequency ui, wavenumber fc, and polarization a. The 
Dirac field operator "L is a linear combination of the an- 
nihilation operators for particles CE>o.k^ E>o,k and the 
creation operators for holes c^, <0 fe ,$ e> <o,k' ■ 

If we now consider the transition matrix elements 
(out I JTint I in) with an initial photon |in) ~ a^, k a\^) an< ^ 
a final particle-hole pair (out) = c ! E>0 k c^ E , <Q k , |0), we 
get to first order in perturbation theory 

1 
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where A° encodes the polarization of the photon. As 
usual, the t- integral gives 6(u> — E + E'), i.e., energy 
conservation. Since the wavelength of the photons un- 
der consideration (in the optical or infra-red regime) is 
much larger than the typical length scales of the elec- 
tronic modes in graphene, we may neglect the photon 
wavenumber k. Therefore, the y-integral yields 8{k — k'), 
i.e., we excite particle- hole pairs with the same wavenum- 
ber k = k' . The remaining a;- integral reads 
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j=E-E', ktzO,<j 
E,k:E',k'=k 



(11) 



Let us first assume AZ = const and consider the tran- 
sition between modes of the same E (i.e., E = —E'), 
such as the upper-most negative mode (for a given k) 
and the lower-most positive mode, cf. Fig. [2] In this 
case, we may use the aforementioned particle-hole sym- 
metry ^ —e k — <jZ ^E,k an d simplify the integrand via 
^emI^^E'M = ^emI^^^em- Inserting 7 1 = ia v and 
7 = — io x and using the properties of the Pauli ma- 
trices, we see that the matrix element for the photon 
polarization in rr-dircction A x yields the same expres- 
sion as in the current J v , cf. Eq. (|5|), and vice versa. 
Consequently, the matrix elements (jllj) vanish for the 
photon polarization in y-dircction, but yield a non-zero 
contribution for the photon polarization in x-direction, at 
least if k is large enough [cf. the discussion after Eq. ([7}]. 
Moreover, the modes with large currents J y and thus 
large group velocities dE/dk do also have large matrix 
elements, which enhances the magneto-thermoelectric or 
magneto-photoelectric effect we are interested in. 

Pseudo-parity Further selection rules arise if we as- 
sume reflection symmetry B(—x) = ~B{x) and thus 
A{— x) = A(x) which yields the additional symmetry 



±iil)2 ,k (x) 



iP E ,k ip2' k (x) ■ 



(12) 



where we call Vem — ±1 the pseudo-parity of this mode. 
Recalling the particle-hole symmetry ^-em = o~ z ^E,k, 
we find V-em = —Vem- The pseudo-parity of a given 
mode can be determined easily for large and positive 



k, where we have i^' ~ vpkipf' k / E from Eq. ([3]). 
Since the wave-function ipf' k (x) of the lowest positive 
mode (for k — > 00) corresponds to the ground state 
of a harmonic oscillator, it is Gaussian and symmet- 
ric ipf' k {— x) — %pf' (x). Hence this mode has an 
even pseudo-parity Vem = +1- The wave-function 
ipf' k (x) of the next mode corresponds to the first excited 
state of a harmonic oscillator and thus is anti-symmetric 
■tpf' k (— x) = — i/jf'(x), which gives an odd pseudo-parity 
Vem = — 1 an d so on. Together with the above result 
V-em = —Vem we that, for a fixed k, the pseudo- 
parity of the modes alternates if we go up and down in 
energy. Assuming that the modes deform continuously if 
k changes [i.e., that A(x) is sufficiently well-behaved], we 
may deduce an alternating pseudo-parity for all k. 



Now, the integrand in the matrix elements (fTTj) be- 
haves as i/j^' (x)tp2 ' k ( x ) i ^ ' k ^ x )'4'i ' k ( x ) f° r the two 
photon polarizations. Inserting Eq. (fT2j) and integrat- 
ing over x, we see that the matrix elements (|11[) be- 
tween modes of the same pseudo-parity vanish for photon 
polarizations in x-direction whereas the transition be- 
tween modes of opposite pseudo-parity is forbidden for 
the other polarization. 

Yet another set of selection rules can be obtained in the 
asymptotic regimes. For k — > 00 we only get transitions 
between modes of opposite energies (due to the orthogo- 
nality of the harmonic oscillator eigen- functions). In the 
opposite limit k — > —00, we recover the well-known (lij 
properties of the Landau levels £7£ = ±vpy/2Bon with 
n£M where we only get transitions for n — > n ± 1. 

Polarization dependence So far, we have discussed 
the case A° = const in Eq. (fTTj) . This is certainly a good 
approximation if the polarization of the incident photon 
points in y direction, i.e., is aligned with the symmetry of 
our set-up. However, for the other (x) polarization, AZ_ in 
Eq. (TiT]) should be replaced by the local projection of the 
photon wave function A a ^ onto the graphene plane, i.e., 
become x-dependent A*(x). The profile of AZ_(x) then 
depends on the incidence angle of the photon. If the 
photon is incident from top, i.e., propagates parallel to 
the external magnetic field k\\B, the two graphene sheets 
(top and bottom) have opposite projections. Thus A a ^{x) 
is anti-symmetric A a ^{— x) = — A°(x) and the above se- 
lection rules are reversed. If the photon propagates per- 
pendicularly through the fold (k _L B), we get a symmet- 
ric projection function A°(— x) = A a ^{x) which vanishes 
far away from the folding region (i.e., for large \x\). In 
this case, the above selection rules do still apply, but the 
matrix elements might be reduced a bit. 

Example profile In order to visualize the behavior of 
the modes by means of a concrete example, let us consider 
a magnetic field of the following form 



B(x) = -Botanh(acc) 



(13) 



where 1/a measures the width of the fold. For a — > 00, 
we get a step function B(x) = i?osign(a;) with the vector 
potential A(x) = -BoM, cf. In this limit, the mode 
equation ([3]) can be solved exactly (piecewise) in terms 
of parabolic cylinder functions, cf. [181 ] . Incidentally, the 
spectrum for such a step function Bjx) = i?osign(cc) can 
also arise for some edge states fl2, 13|. 

However, such a step function can only be a good ap- 
proximation if k is not too large and if the curvature 
radius of the graphene fold is much smaller than the typ- 
ical magnetic length scale £3 = l/v^B- For one Tesla, 
we get £b ~ 26 nm while the radius of curvature cannot 
be too small since it should be much larger than the lat- 
tice spacing ps 0.25 nm. Thus, let us consider a finite a 
and take a = as an example. The spectrum can 

then be obtained numerically and is given in Fig. [5J The 
spectra for other values of a are qualitatively similar. As 



4 



demonstrated above, the two lowest modes are monoton- 
ically increasing/decreasing, whereas the higher modes 
can have dE/dk = at some small k- values. For large 
| A; |, we recover the asymptotics discussed above. 
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FIG. 2. Dispersion relation of the lowest bands with k = Mb 
and E = evf\/2Bo and sketch of the photo-absorption. 

Conclusions Via the effective Dirac equation ([1}, we 
studied the low-energy behavior of electronic excitations 
in graphene under the influence of a transversal magnetic 
field B{x) with the asymptotics B(x — > ±oo) = ±Bq. 
Such a field profile B(x) arises within a folded graphene 
sheet in a constant magnetic field, for example, see Fig.[T] 
Based on general arguments, we find a discrete set of 
modes (see also [20| ) which are localized near the fold 
(i.e., the zero of the magnetic field) and propagate along 
it with a significant fraction of the Fermi velocity. 

Due to particle-hole symmetry, the dispersion relations 
E(k) of these modes (cf. Fig. [5J are symmetric around 
the E = axis, but do never cross it. Thus, these modes 
have a finite energy gap (for each k) with the character- 
istic energy scale being set by the (pseudo-relativistic) 
Landau level energy El = vpy/2Bo. For a magnetic field 
of one Tesla, we have El ~ 36 meV, which corresponds 
to 400 Kelvin. The group velocity dE/dk is related to the 
current J v and we find that particles and holes move in 
opposite directions. Apart from some minor exceptions, 
all particles move to the right and all holes move to the 
left - i.e., we get a nearly perfect charge separation. In 
view of this pre-determined direction, the finite energy 
gap, and the fact that these localized modes are qualita- 
tively independent of the shape of A(x), we expect that 
they are quite robust against perturbations. In addition, 
we consider the propagation within a (curved) graphene 
sheet, i.e., far away from any edges with defects etc. 

Finally, we study the excitation of particle-hole pairs 
in these modes via incident infra-red or optical photons, 
i.e., the magneto-thermoelectric (Nernst-Ettingshausen) 
or magneto-photoelectric effect. The matrix elements 
display a distinct dependence on the polarization and the 



incidence angle of the photons, which should enable us to 
distinguish this effect from other phenomena experimen- 
tally. Furthermore, we find that those modes with com- 
parably large group velocities (i.e., large currents) tend 
to have large matrix elements (at least for low-energy 
transitions) and thus are more strongly coupled to the 
incident photons (i.e., "nature favors our goal"). 

Outlook: electric field If we apply an additional elec- 
tric field perpendicular to the fold and the magnetic field, 
we get an electrostatic potential $(x) = f3vpA(x) with 
some constant f3. If we have \/3\ < 1 (i.e., if the elec- 
tric field is sub-critical), we may transform $ away by 
an effective Lorentz boost in y-direction with a veloc- 
ity ^hnnst, = ftv-F where v-p plays the role of the speed of 
light [19( . In the Lorentz boosted frame, we get the same 
modes as discussed above, but with a reduced magnetic 
field B' = Bq y/l — /3 2 . Since this field enters the charac- 
teristic energy scale via vf\/2Bq, the dispersion relation 
after transforming back to laboratory coordinates reads 

E^E' = £(l-/3 2 ) 3 / 4 -/h; F /3, (14) 

i.e., the spectrum in Fig. [2] is compressed and tilted. 
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